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ABSTRACT 

In this paper we present detailed analysis of the transient X-ray source 2XMMi 
J003833.3+402133 detected by XMM-Newton in January 2008 during a survey of M31. This 
source has previously been identified as a black hole X-ray binary in M31, but here we argue 
that the X-ray spectra and timing data are inconsistent with this conclusion. We instead argue 
that 2XMMi J003833. 3+402133 may be a new addition to the rare class of magnetars. The 
X-ray spectrum is well fitted by either a steep power law plus a blackbody model or a double 
blackbody model. Prior observations with XMM-Newton, Chandra, Swift and ROSAT spanning 
1991 to 2007, as well as an additional Swift observation in 2011, all failed to detect this source. 
No counterpart was detected in deep optical imaging with the Canada France Hawaii Telescope 
down to a 3cr lower limit of g = 26.5 mag. The transient behaviour, X-ray spectrum, and lack 
of an optical counterpart are all consistent with a magnetar. The derived luminosity and black 
body emitting radius at the distance of M31 argue against an extragalactic location, implying 
that it is located within the Milky Way but 22° out of the plane. The high Galactic latitude 
could be explained if 2XMMi J003833. 3+402133 were an old magnetar, or if its progenitor was 
a runaway star that traveled away from the plane prior to going supernova. 

Subject headings: Stars: Neutron, X-rays: individual (2XMMi J003833. 3+402133), X-rays: stars 



1. Introduction 

Magnetars are a rare class of neutron stars 
which are characterised by a very strong sur- 
face magnetic field, typically 10 14 - 10 15 G 
dDuncan fe Thompson 1992t Woods fe Thompson! 



gamma rays. While originally considered to be un 
related, it is now accepted tha t AXPs and SGRs 
are the same type of object (iKouyeliotou et al 



2006), and relatively slow spin rate (P ~2 - 12 s) 
when compared to the majority of observed pul- 
sars. They also differ from normal pulsars in that 
they are not powered by their spin-down energy 
losses but by t he decay of energy stored i n their 
magnetic field (jDuncan fe Thompsonlll992l ). 

Magnetars are historically divided into two 
classes: anomalous X-ray pulsars (AXPs) which 
are characterised by slow pulsations and rapid 
spin-down, and soft gamma repeaters (SGRs) 
which exhibit sporadic bursts of hard X-ray/soft 



1998[ iGavriil. Kaspi. fc Woods! 120021 iKaspi et al 
2003). There are currently 17 confirmed mag- 
netars (8 SGRs and 9 AXPs) and a further 5 
candidates (2 SGRs and 3 AXPsfl- Due to 
this small sample size it is difficult to deter- 
mine the average lifetime, spatial distribution, 
duty cycle and luminosity of the population as 
a whole. However, of the known magnetars the 
majority lie within the Galactic plane, consistent 
with young (<10,000 yrs) neutron stars with low 



1 Prom McGill's Online Catalogue available at 
http: //www. physics .mcgill . ca/~ pulsar /magnetar /main, html 
with the addition of the recently confirmed SGR 1822-1606 
jGogiis et alj|201lf) . 



spati al velocities (< 500km s 1 ; iGaensler et al 
20011 ). The exceptions are the t wo magnetars 



locat e d in the Magellanic Clouds (jMazets et al 
19791 : iMcGarrv etail 120051) . While there have 



been reports of the detection of SG Rs in M31 



Ofck et all 2008: Ma zets et all 120081 ) and M81 



Frederiks et al.l 20071 ). these sources have been 



detected only once and their SGR nature has thus 
yet to be verified. 

In this paper we present an in-depth analy- 
sis of multi-wavelength data on the X-ray source 
2XMMi J003833.3+402133 (hereafter, referred to 
as XMM J0038+40), which was obse rved during 



a sury ey of the nearby galaxy M31. IStiele et al 
( 201l[ ) identified this object (source 57 in their 
sample) as a probable black hole low mass X-ray 
binary (LMXB) located within M31. Here we ar- 
gue that the X-ray spectral and timing properties 
of this source are not consistent with a black hole 
LMXB, and instead propose that XMM J0038+40 
may be a new addition to the rare class of Galactic 
magnetars. If this conclusion is correct, the high 
Galactic latitude of XMM J0038+40 would make 
it unique as the other Galactic magnetars are all 
located within the plane. 

2. Data Reduction, Analysis &: Results 
2.1. X-ray Data 

The field of XMM J0038+40 was observed by 
XMM-Newton on 2008 January 2 for 45.5 ks by 
the three EPIC cameras (pn, MOS1 and MOS2) 



as part of a 


survey 


(Stiele et al. 


2011) 



Source detection performed 
by the XMM-Newton processing pipeline detected 
a bright uncatalogued source at an off-axis angle of 
7.72'. The source position is RA = 00h38'33.32", 
dec = +40°21'33.20" (J2000 ) with a 3a positiona l 
uncertainty radius of 1" (see lWatson et al. ( 20091 ) 
for a detailed description of the source detection 
procedure). No data were available from the two 
Reflection Grating Spectrometers (RGS) as XMM 
J0038+40 fell outside their fields of view. 

The Science Analysis System (SAS) version 
10.0 softwar^l was used to process the Observa- 
tion Data Files (ODFs) usin g the s ame procedure 
as outlined by iFarrell et al] (|2010l) . Background 
light curves extracted from the entire field at ener- 



gies > 10 keV found no evidence for flaring, and so 
no good time interval filtering was applied. We ex- 
tracted source spectra using circular regions cen- 
tred on the source position with radii of 71.5" and 
41.5" for the pn and MOS cameras, respectiveljU. 
Background spectra were extracted from nearby 
source- free regions of the field with areas 1.5 times 
the source extraction region areas. Response and 
ancillary response files were produced for the spec- 
tra and the light curves were corrected using the 
SAS task epiclccorr. The spectra were grouped at 
20 counts per bin to provide sufficient statistics 
for y 2 fitting a nd analysed using XSPEC vl2.6.0q 
(lArnaudl[l996h . ^)ata below 0.2 keV and above 
7 keV (where the statistics are very poor) were 
ignored for the spectral fitting. 

2.1.1. X-ray Spectral Analysis 



Stiele et al.l (|2011l ) fitted spectra extracted from 
their 2008 XMM-Newton observation with three 
models: absorbed disk blackbody plus power 
law, absorbed disk blackbody, and absorbed 
bremsstrahlung models. They found that the 
disk blackbody plus power law model was the 
best fitting model obtaining % 2 /dof = 174/145, 
although they also claim formally acceptable 
fits using the two other models (with x 2 /dof = 
270/147 and 209/147 for the disk blackbody and 
bremsstrahlung models, respectively). 

We fitted the spectra we extracted from the 
2008 XMM-Newton observation with the same 
models and obtained similar results for the disk 
blackbody plus power law and the simple disk 
blackbody models, although we argue that due to 
the high quality of the data the fit with the sim- 
ple disk blackbody model is unacceptable with a 
X 2 /dof = 812/714. In contrast, our fit with the 
absorbed bremsstrahlung model obtained a much 
better fit (x 2 /dof = 764/714) with a significantly 
lower temperature of kT = 0.98 ± 0. 02 keV than 
the kT = 1.91 ± 0.07 keV obtained bv lStiele et ail 



(|201ll ). The cause of the di fferences between o ur 
fitting and that reported bv lStiele et al.l ( 2011 ) is 
unclear. We also fitted the spectra with a more 
physical model representing emission from an ac- 
creting black hole, i.e. an absorbed disk blackbody 



2 http: //xmm. esa. int/sas/ 



3 The source extraction radii were chosen so as to maximise 
the signal-to-noise ratio of the data products (V. Braito 
2011, private communication). 
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(Ti n — 0.30 keV) plus thermal Comptonisation 
model (compTT in XSPEC). Again, we obtained 
an acceptable fit with x 2 /dof = 748/711 and spec- 
tral parameters of T = 0.31 keV, kT = 0.49 keV, 
and t = 23.79. 

In addit i on to the model fits reported by 



Stiele et al.l (|201lh . we attempted to fit the 2008 



XMM-Newton spectra with power law, black body 
(the BBODYRAD model in XSPEC), and ther- 
mal plasma (the MEKAL model in XSPEC) 
models. In each case photoelectric absorp- 
tion was accounted for using the phabs com- 
ponent in XSPEC and the W ilms abundances 
l| Wilms. Allen, fc McCravll2000l ) . Neither the sim- 
ple power law model nor black body models pro- 
vided an acceptable fit, with x 2 /dof = 1205/719 
and 1165/719, respectively, and significant resid- 
uals appearing below 2 keV. Adding a low tem- 
perature black body component to the power law 
model improved the fit significantly (x 2 /dof = 
806/712), although with a very steep power law 
photon index (see Table [lj . The addition of a sec- 
ond black body component to the simple absorbed 
black body model produced a slightly better fit 
(x 2 /dof = 752/712) and completely smoothed out 
the low energy residuals (see Table Q] and Figure 

Attempts to fit the spectra with an absorbed 
MEKAL model with the abundance parameter 
frozen at Solar values did not provide an accept- 
able fit (x 2 /dof = 11527/719). Allowing the abun- 
dance to vary freely improved the fit significantly 
(x 2 /dof = 769/718), however the abundance value 
fell to zero indicating that no significant line emis- 
sion is present and therefore the model is consis- 
tent with the underlying bremsstrahlung contin- 
uum model. In summary, we obtained acceptable 
fits with the double blackbody, blackbody plus 
power law, disk blackbody plus power law, and 
bremsstrahlung models but are able to rule out 
the simple black body, simple power law, and ther- 
mal plasma (with non-zero elemental abundances) 
models. 

2.1.2. X-ray Timing Analysis 

To search for pulsations in the XMM-Newton 
data, we extracted counts in the energy range 0.2 
- 12 keV using the same extraction regions as used 
for the spectra. This resulted in 8078, 7192 and 
23740 counts for the MOS1, MOS2 and pn cam- 



eras, respectively. In each case, we corrected the 
photon arrival times to the Solar system barycen- 
ter. To search for high frequency periodic variabil- 
ity we performed a Z\ (Ra yleigh) test to sear ch for 
sinusoidal pulsations (e.g. iLeahv et al.lll983l) . For 
the MOS1 and MOS2 data, we searched for pulse 
periods in the range 6.4 to 4500 s, while for pn data 
we considered periods between 146 ms and 4500 s. 
We saw no significant power at any period, either 
in each data set separately or in a period search 
of all three data sets combined. The correspond- 
ing 5cr upper limit on the pulsed fraction is 7% in 
the period range 146 ms to 6.4 s, and 5.5% in the 
period range 6.4 to 4500 s. 

We also extracted power density spectra (PDS) 
over the reliable energy bandpass (0.3-6 keV) 
binned on 10 s. We find that, with linear bin- 
ning, the variability is consistent with the sta- 
tistical white noise within \a using the most 
accurate Bayes ian least-squares fitting routines 
(|Vaughanll2010h . As the shape of the PDS is en- 
ergy dependent and we could conceivably have a 
stable, bright component diluting the variabil- 
ity when averaged across the entire bandpass, 
we also extracted the ener gy dependent frac- 
tional excess variance (see e.g. lEdelson et al.ll2002t 
IVaughan et al.l 120031 ) . We did this initially for 
the bands which correspond to the absorbed disk 
blackbody plus thermal Comptonisation spectral 
model (0.3 - 0.8 and 0.8 - 6keV; see Figure W) 
with binning on 25, 50, 100, 200 and 400 s to test 
for variability typically seen from black hole X-ray 
binaries. We could not constrain any variability 
over any timescales in either energy band above 
the white noise level. As this relies on our spectral 
deconvolution being an accurate description of the 
data we also tested energy bands within these two 
bands (specifically 0.3 - 0.5, 0.5 - 1, 1 - 2 and 
2-6 keV) on 400 s binning (to ensure adequate 
statistics in these smaller energy bands) and find 
no constrained variance above the statistical white 
noise. 

2.1.3. Additional X-ray Observations 

M31 has been observed numerous times pre- 
viously in X-rays, and thus archival data from 
ROSAT, XMM-Newton, Chandra, and Swift were 
also analysed. In addition, following our discovery 
of XMM J0038+40 we obtained Target of Oppor- 
tunity (ToO) observations of the field with the 
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Fig. 1.— XMM-Newton EPIC spectra (black 
= pn, red = MOS1, green = MOS2) of XMM 
J0038+40 from the 2008 detection fitted with the 
Phabs*(BB+BB) model. The bottom panel shows 
the fit residuals. 
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Fig. 2. — Spectral model components for the 
absorbed disk blackbody (red line) plus thermal 
Comptonisation (blue line) fit. 



Swift X-ray Telescope (XRT) between 04-15 
February 2011. For the additional XMM-Newton 
data, we extracted images from all three EPIC 
cameras after filtering out background flares. 
Pre-processed images from the archival Chandra, 
ROSAT, and Swift data were downloaded from 
NASA's HEASARC archive... The 2011 Swift ToO 
observation was processe d using the online XRT 
data processing facilitjH ([Evans et al.ll2009t ). 

XMM J0038+40 was not detected in any of 
these data (Figure[3]). Upper limits were estimated 
using the count rates measured within regions cen- 
tred on the position of XMM J0038+40 in each 
of the images, subtracted from background count 
rates measured in nearby source-free regions. The 
WebPIMMS online flux conversion toojj and the 
double blackbody spectral model obtained from 
fitting the EPIC spectra from the 2008 XMM- 
Newton observation were used to convert the 3<r 
count rate upper limits into flux limits. The flux 
limits derived in this manner are all well below the 
flux observed in the 2008 XMM-Newton observa- 
tion indicating that XMM J0038+40 is a transient 
source with variability by a factor > 450. Specific 
details of the observations from which data was 
analysed including the upper limits of the non- 
detections are given in the online Appendix (Table 

2.2. UV, Optical and Near-Infrared Data 

XMM J0038+40 fell within the field of view of 
the Optical Monitor (OM) during all three XMM- 
Newton observations, which observed the field in 
imaging mode with the uvwl and uvm2 UV fil- 
ters. As the 2008 OM images are the only available 
UV data that were taken at the time when XMM 
J0038+40 was X-ray bright, we focused on these 
data for our analysis. The total exposure times 
during this observation were 1461 s and 6317 s for 
the for uvwl and uvm2 filters, respectively. No 
source was detected within the 3ct X-ray error cir- 
cle of XMM J0038+40 in any of the images. We 
determined 3<r count rate upper limits at the po- 
sition of XMM J0038+40 using the XMM-Newton 
pipeline processed images and the same method as 
that used for the X-ray non-detections, and con- 



4 http: //heasarc .gsf c .nasa. gov/cgi-bin/W3Browse/w3browse .pi 

5 http: //www. swif t . ac.uk/user_objects/ 

e http: //heasarc .gsf c .nasa. gov/Tools/w3pimms .html 



4 



Fig. 3.— X-ray images of the field of XMM J0038+40 taken with XMM-Newton (panels 1 to 3 from the 
left) and Swift. The white circle indicates the position of XMM J0038+40 with a radius of 30". The XMM- 
Newton images span the energy range of 0.2 - 12 keV, while the Swift image spans 0.3 - 10 keV. A Gaussian 
smoothing function with a kernel radius of 3 pixels was applied to all images. 



verted the count rates into Vega magnitudes using 
the filter zero points for the 

Deep g- and i-band images of the field of XMM 
J0038+40 were taken with the Canada France 
Hawaii Telescope (CFHT) on 2008 October 3 as 
part of the Pan- Andromeda Archaeo logical Survey 
fPAndAS: [McConnachie et alJl2009h . The CFHT 
observations and data r e ductio n procedures are 
described by llbata et al. ( 2007 ). No counterpart 
was detected in either band image within the X- 
ray error circle of XMM J0038+40 (see Figure gj}, 
so we determined magnitude lower limits in the 
same manner as described above. There is how- 
ever a faint point source just outside the error 
circle at RA = 00h38'33.3", dec = +40°21'31.9" 
with g = 24.3 ± 0.1 mag and i = 22.11 ± 0.04 
mag which could possibly be the counterpart. No 
source was present within or even nearby the X- 
ray error circle in the 2MASS J, H, and K s images. 
In order to estimate conservative magnitude lower 
limits in the 2MASS near-infrared (NIR) data, we 
identified the faintest source in each band within 
30' of XMM J0038+40 and adopted these mag- 
nitudes as the upper limits. All the lower limits 
derived from the UV, optical and NIR data are 
given in Table [2] 

2.3. Radio Data 

The region surrounding XMM J0038+40 has 
been imaged at various radio frequencies as part 
of surveys of M31. Using these data, we find no 
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Fig. 4. — CFHT g-band image with the X-ray po- 
sition of XMM J0038+40 indicated by the white 
circle, with a radius equal to the 3<r positional er- 
ror of 1". 



Table 2: UV, Optical and NIR Lower Limits 



7 http : //xmm. esa. int/ext ernal/xmm_user_support /document at ion/uhb_^ . b/node I b .htm± 
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H 
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radio source at the position of XMM J0038+40 
at frequencies of 0.3, 0.6 or 1.4 GHz, with 5cr up- 
per limits on the flux of a ra dio counterpart of 5 , 
4 and 1.5 mJy, respectively (IBvstedt et al.l Il984 ; 



Walterbos et all 11985: Gcl fand et all 120041 ). Less 



sensitive observations at 0.3 and 1.4 GHz at other 
epochs also show no radio source at this p osition 
(|Rengelink et al.lll997t ICondon et alJll998t ). 



3. Discussion &: Conclusions 

The X-ray spectrum, transient behaviour, 
and lack of an optical and infrared counterpart 
place strong constraints on the nature of XMM 
J0038+40. Non-transient objects such as isolated 
cooling neutron stars and rotation powered pul- 
sars with moderate magnetic field strengths can 
immediately be ruled out, leaving behind objects 
that are known to exhibit large scale variability. 

Although the CFHT optical imaging was not 
taken at the same time as XMM J0038+40 was 
detected in X-rays, we would still expect to detect 
an optical counterpart consistent with the X-ray 
source position if XMM J0038+40 were a star, a 
cataclysmic variable (CV) or an Active Galactic 
Nucleus (AGN). The X-ray detection and deep g- 
band and z-band limits indicate X-ray to optical 
flux ratios of F x /F g > 8200 and F x /F f > 5800 for 
XMM J0038+40. Even if the source just outside 
the XMM-Newton error circle were the counter- 
part the flux ratios would be Fx/F g ~ 1300 and 
Fx /Fj ~ 340. These are far in excess of the ratios 
observed for AGN and stars, which typically have 
X-ray to optical flux rat ios of Fy/Fp < 10 and 
< 0.01, respectively (e.g. iMainieri et al.lliooH ). It 
is also inconsis tent with CVs which typically have 
Fjf/Fo < 0.1 feuulkers et al.ll2006l) . 

High levels of dust extinction could explain the 
lack of an optical counterpart!!, however the posi- 
tion of XMM J0038+40 -22° out of the Galactic 
plane puts it in a region o f low dust extinctio n 
with E(B-V) = 0.062 mag (jSchlegel et al.lll998l) . 
The weighted average Galactic neutral hydrogen 
absorption in the direction of XMM J0038+40 as 
measured by the Leiden/Argentine/Bo nn (LAB) 
surve y of Galactic HI is 7.6 x 10 20 cm -2 (jKalberla et al 
20051 ). consistent with the N# values determined 
through the X-ray spectral fitting with some ad- 



'But not the lack of a NIR counterpart in the 2MASS im- 
ages. 



ditional intrinsic absorption. If XMM J0038+40 
is located within or behind M31 there may be 
an additional contribution from the galaxy disk. 
However, the low N# values from the X-ray spec- 
tral fitting indicate that M31 is unlikely to be a 
major contributor to the extinction. We thus con- 
clude that the high flux ratio is not due to dust 
reddening. 

The lack of a point-like radio counterpart above 
1.5 mJy is not very constraining. For example, us- 
ing the black hole fundamental plane relationship 
between X-ray luminosity, radio luminosity, and 
black hole mass, we wouldn't expect to be able to 
detect radio emission from an AGN unless it was 
closer than 1.5 M pc (assuming a bla ck hole mass 
of 10 6 M Q ; e.g. iKording et "all 120061) . At such a 
distance the galaxy itself should be easily resolved 
in the optical/NIR images. Likewise, radio emis- 
sion from the jet of a 20 M black hole X-ray 
binary would not be detectable unless the system 
was closer than 20 pc, a distance at which emis- 
sion from the disc or donor star would surely be 
detectable in our deep optical images. 

The EPIC X-ray spectra are consistent with 
black body plus steep power law, double black 
body, simple disk blackbody, disk blackbody plus 
power law, or thermal bremmsrahlung models (all 
with the addition of a photoelectric absorption 
component of low column density, consistent with 
the Galactic latitude of the source) . The temper- 
ature of the bremmstrahlung model was 0.98 keV, 
much lower than the lowest temperature measured 
for a CV (kT = 1.9 keV from thermal plasma fits 
to XMM-Ne wton spectra of th e quiescent dwarf 
nova GW Lik lHilton et al.ll2007l ). In addition, line 
emission becomes more significant relative to the 
underlying continuum spectrum for CVs with low 
temperature spectra, and so if XMM J0038+40 
were a CV we would expect that the spectrum 
would be consistent with a thermal plasma model 
with non-zero elemental abundances. The spec- 
trum of XMM J0038+40 is not consistent with 
such a model, adding further weight to the argu- 
ment against it being a CV. 

Spectra of accreting black hole and neutron 
star X-ray binaries are typically modelled by a 
power law (representing inverse Compton emis- 
sion), sometimes with the addition of a soft ther- 
mal component attributed either to emission from 
the surface (in the case of a neutron star) or the 
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accretion disk (e .g. iRemillard fc McClintocki r2006: 
OrlandinillioOffl . Observationally it can be diffi- 
cnlt to discriminate between black hole and neu- 
tron stars in X-ray binaries unless the compact ob- 
ject undergoes behaviour such as thermonuclear 
bursts or spin period pulsations that are defini- 
tively identified with neutron stars. However, the 
inverse Compton tail in the X-ray spectra of black 
hole X-ray binaries is typically much steeper than 
that seen in neutron star X-ray binaries, allow- 
ing them to be separated via their X-ray spec- 
tral hardness (Farrell et al. 2011, in preparation). 
When the EPIC spectra of XMM J0038+40 were 
fitted with the blackbody plus power law model, 
the slope of the latter component was very steep 
with a photon index of 3.7, more consistent with 
black hole binaries than neutron stars (although 
significantly steeper than typically seen even for 
black hole binaries in the high luminosity spectral 
states). Replacing the blackbody component with 
a multi-color disk blackbody model resulted in a 
flatter power law slope, although still steeper than 
seen for neutron star X-ray binaries. 

The temperatures derived for the disk black- 
body components when fitting the EPIC spectra 
either with the simple disk blackbody model or 
the disk blackbody plus power law model are sig- 
nificantly less than that generally seen in black 
hole binary spectra when the disk is inferred 
to be at the inner sta ble circular orbit (ISCO ; 
see t he reviews 
20061: 



of 



Done et al.l [2007). 



McClintock fe Remillard 



If the disk is trun- 
mass accretion rates (see 



cated at t he low est 

Maccarond 120031 ) then the temperature of this 



component should be lower and the spectrum 
shoul d be appended by a strong, har d tail of emis- 
sion (jMcClintock fc Remillardl 120061 ) . This tail is 
always seen with large amounts of variability on 
all timescales (although this variability itse lf may 



originate in the disk: see lUttlev et al.ll201ll ). The 



shape of the X-ray spectral data alone do not sup- 
port an identification with such an accretion state 
and we find concurrent evidence from inspecting 
the variability. This strongly suggests that the 
observed emission is not associated with accretion 
onto a black hole at low fractions of the Edding- 
ton luminosity, and yet the disk temperature is 
too low for accretion onto a black hole at a high 
Eddington ratio. 

We note that placing the source at ~2 kpc 



results in a luminosity consistent with that of a 
black ho le X-ray binary in q uiescence (~ 10 33 erg 
s _1 ; e.g. lBradlev et al.1120071 ). where the spectrum 
should be dominated by an advection dominate d 
accretion flow (ADAF; see Naravan et al.l [l997T) . 
The best studied of such systems is V404 Cyg (e.g. 
Bradley et ap2007h which has a soft Comptonised 



component (with T >2) although this may have 
a contribution from the disk (now truncated at 
large radii). This matches the 2008 XMM- Newton 
EPIC spectral data for XMM J0038+40 extremely 
well. However, the ADAF is still ex pected to be 
signifi cantl y variable. Fo r exam ple, iHvnes et al.l 
(l2004h and lBradlev et al.1 <|2007h found V404 Cyg 
in quiescence to be highly variable in light curves 
binned at 400 s and 200, respectively. As the 
ADAF totally dominates the emission in the quies- 
cent state the lack of variability over the full band- 
pass argues strongly against a quiescent black hole 
binary. 

In addition, a distance of ^2 kpc would imply a 
high state flux of ~ 10~ 6 erg cm~ 2 s _1 (assuming 
a 10 Mq black hole emitting isotropically at 80% 
of the Eddington limit) which would make XMM 
J0038+40 one of the brightest X-ray sources in the 
sky during the high state. The lack of a detection 
in any archival X-ray data that we examined or 
in the regular all sky surveys performed with the 
RXTE satellite is thus difficult to conceive in this 
scenario, and implies that XMM J0038+40 was in 
the high state as opposed to the quiescent state 
during the 2008 XMM-Newton observation. We 
therefore conclude that XMM J0038+40 is not a 
black hole X-ray binary. 

The X-ray spectra of magnetars have been 
historically fitted with a two component model 
comprised of a steep power law p lus a black- 
body (e.g. Marsden fe White! 120011) . However, 



Halpern fc Gotthelfl (|2005) argue that this model 



is not physically justified and favour the fitting of 
a double blackbody model. They emphasize this is 
more physically feasible as the pulse modulation 
of the magnetar can be understood in terms of 
anisotropic heat conduction and radiative transfer 
in the strong magnetic field. The EPIC spectra 
of XMM J0038+40 are consistent with both the 
blackbody plus steep power law and the double 
blackbody models, with the spectral parameters 
falling well with in the range spanned by SGRs 
and AXPs (e.g. iMarsden fc White! l200lh . The 
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absence of slow pulsations in the light curve of 
XMM J0038+40 (one of the defining character- 
istics of magnetars) does not weigh against the 
magnetar hypothesis, as pulsed fractions below 
the upper limit of 7% derived for XMM J0038+40 
have been observed from other magnetars (e.g. 
the AXP 4U 0142+614 has be en observed with a 
pulsed fraction as low as ~5%: lRea et~aDl2n07h .. 

Magnetars typically have very faint or no coun- 
terparts in NIR and optical wavelengths (e.g . 
iHulleman et "all 1200(1 iGelfand fc Gaenslerll2007i 
consistent with the deep limits we have obtained 
for XMM J0038+40. In particular, using the un- 
absorbed X-ray flux and the 2MASS K s band 
magnitude limit, the X-ray to NIR flux ratio for 
XMM J0038+40 is F x /F Ks > 37, which is consis- 
tent with the flux rati os of the known magnetar s 
(see e.g. Figure 5 of IGelfand fc Gaenslerl 120071) . 
The X-ray variability by a factor of > 450 observed 
from XMM J0038+40 is of the same order as 
seen from other magnetars (e.g. XTE J1810-197, 
CXOU J164710.2-455216 and AX J1845.0- 0258; 
see e.g. Tarn et all 12006 ; lAlbano et al. 2010l and 
references therein) but slightly higher than the 
most extreme example (AX J1 845. 0-0258 with a 
variability factor of 260 - 430; iTam et al.ll2006l) . 
Thus, the X-ray spectrum, transient behaviour, 
and high Fx /Fq ratio are however all consistent 
with the known sample of magnetars, leading us 
to conclude that XMM J0038+40 is most likely a 
new addition to this rare class of object. 

Apart from the two magnetars located in the 
Magellanic clouds all the other members of this 
class have been found within the plane of our own 
Galaxy. In contrast, XMM J0038+40 is coinci- 
dent with the outer edges of the nearby galaxy 
M31. The 0.2 - 10 keV X-ray luminosity of XMM 
J0038+40 at the distance of M31 (D ~ 0.8 Mpc) 
is ~ 5 x 10 37 erg s _1 , about 2 orders of magni- 
tude higher than typi cally seen from magn etars 
following an outburst (|Rea fc Espositoll201lf) . In 
addition, if XMM J0038+40 is located in M31, 
the radius of the black body emitting radii de- 
rived from the normalization of the BBODYRAD 
components are 554 km and 126 km for the double 
black body model and 162 km for the black body 
plus power law model. These radii are far greater 
than the maximum radius of a neutron star. When 
taken in conjunction with the unusually high lu- 
minosity, we conclude that XMM J0038+40 is un- 



likely to be a magnetar located within M31, and 
is therefore most likely a new Galactic magnetar. 

If XMM J0038+40 is < 10,000 yr old with a spa- 
tial velocity < 500 km s _1 , assuming it was born 
in a supernova explo sion from a massive st ar in the 
Galactic plane (e.g. iGaensler et al.l 120011 ). then it 
could only have travelled ~5 pc directly out of the 
plane. The Galactic latitude of b ~22° therefore 
implies a line-of-sight distance of ^13 pc, making 
XMM J0038+40 an improbably nearby magnetar. 
Furthermore, the derived luminosity at such a dis- 
tance is only ~ 5 x 10 27 erg s _1 , significantly lower 
than typically observed from magnetars. If the 
luminosity of XMM J0038+40 is between 10 32 - 
10 35 erg s _1 like the other magnetars, this implies 
a distance of between ~0.7 - 24 kpc. This then 
implies either an implausibly high spatial velocity 
> 10,000 km s" 1 (for an age < 10,000 yr) or an 
age > 3 x 10 6 yr (for a spatial velocity < 500 km 
s- 1 ). 

Such an age is inconsistent with that derived 
for magnetars. However, SGR-like behaviour has 
been observed from a source (SGR 0418+5729) lo- 
cated at a moderate Galactic latitude (b = 5.1°). 
Pulsations with a period of ~9.1 s have been de- 
tected from SGR 0418+5729 but with a very small 
period derivative of P < 6 x 10~ 15 s s _1 , im- 
plying a surface dipolar magnetic field strength 
of < 7.5 x 10 12 G and a chara cteristic age of 

" 20101) . Thus, if the 



> 24 x 10 6 yr (|Rea et al 



progenitor of XMM J0038+40 was a massive star 
that went supernova within the Galactic plane, it 
is possible that it is an old magnetar where the 
magnetic field strength has decayed to ~ 10 12 G. 

Alternatively, the high Galactic latitude of 
XMM J0038+40 could be explained by the pro- 
genitor being a massive runaway star that origi- 
nated in a binary system that became unbound, 
and then travelled outside the plane of the Galaxy 
prior to going supernova. Runaway O and B-type 
stars have been observed in our Galaxy with space 



veloc ities up to ^100 km s ijHoogerwerf et al 



20011 ). Assuming a similar space velocity, the 



progenitor of XMM J0038+40 would have taken 
- 6.8 x 10 6 - 2.3 x 10 8 yr to travel ~0.7 - 24 
kpc. These timescales are perfectly consistent 
with the nuclear burning lifetimes of massive stars 
with masses between ~8 - 30 M^ with lifetimes 
betw een 
19891 ). 



7 - 40 x 10 6 yr (jMaeder fc Mevnet 
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In order to discriminate between these two sce- 
narios we first need to confirm the magnetar na- 
ture of XMM J0038+40 through the detection of 
pulsations. The measurement of the spin deriva- 
tive would then allow us to derive both the surface 
dipole magnetic field strength and the character- 
istic age of the neutron star. Further deep X-ray 
observations triggered based on regular monitor- 
ing of XMM J0038+40 are thus required. 
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Table 1: Best fit spectral parameters to the 2008 XMM- Newton observation of XMM J0038+40. 



Parameter 


Phabs*(BB+BB) a 


Phabs*(BB+Pow) a 


Units 


N ff 


0.08±0.01 


S5 +U U2 
u - 00 -o.oi 


x 10 22 cm" 2 


kTr 


0.207±0.005 


0.37±0.01 


keV 


Norm! 

kT 2 /r b 

Norm2 


n 44+ 01 
2 5 +0 - 7 


4 1+ 1 - 

q 7+O.2 
°-'-0.1 

7.8±g;l x 10- 4 


kcV/ ••• 


Flux abs c 


i-^y_o.oi 


1.5±0.02 


x 10~ 12 erg cm -2 s _1 




1.69±0.02 


i2±| 


x 10~ 12 erg cm -2 s _1 


Fluxi/Flux 2 d 
X 2 /dof 


1.07 

752/712 


0.07 
806/712 





"All errors are quoted at the 90% confidence level. 

''This row contains the kT for the second BB component in the BB+BB model and the photon index of the power law component 
in the BB+Pow model. 
c Absorbed and unabsorbed fluxes are calculated over the energy range 0.2 - 10 keV. 
d Ratio of unabsorbed 0.2 - 10 keV fluxes of model component 1 over model component 2. 
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Details of X-ray Observations 



Table 3 
Log of X-ray Observations. 



Telescope ObsID Date Instrument T ea; p a Flux Upper Limit b 



ROSAT 

± t kj lj /l J. 


600079N00 

KJKJKJKJ 1 i/l > KJ KJ 


1991- 


■07-14 


PSPC 


41.7 


1.7 




10~ 


14 


ROSA T 

J. L KJ kJJi J. 


60oo64noo 

KJ \J \J i / w t: i > ww 


1991- 


-07-15 


PSPC 


48.8 


2.8 


x 


10~ 


14 


ROSA T 


1 41 836noo 

_l_T:_l_L}tJWl , l KJKJ 


1992- 


-01-12 


PSPC 


1.7 


6.7 


x 


10~ 


14 


EOS AT 


1 41 839N00 


1992- 


-01-13 


PSPC 


1.8 


3.6 


x 


10~ 


14 


ROSAT 

J. t KJ kJ Ix. _L 


141837N00 


1992- 


-01-13 


PSPC 


2.1 


1.2 


x 


10~ 


13 


ROSAT 


60031 7Non 

KJKJKJkj -L I liUU 


1992- 


■07-20 


PSPC 


2.6 


5.4 


x 


10~ 


14 


ROSAT 

1 1 KJ kJ jTL J_ 


600301 nod 

UUUOUll> KJKJ 


1992- 


■07-22 


PSPC 


2.7 


3.9 


x 


10~ 


14 


EOS AT 


600341N00 


1992- 


■07-23 


PSPC 


2.7 


9.3 


x 


10~ 


14 


EOS AT 


600360N00 


1992- 


-07-23 


PSPC 


2.9 


1.2 


x 


10~ 


13 


EOS AT 

i ( y y ky ii J- 


600339N00 


1992- 


-07-25 


PSPC 


2.8 


7.9 


x 


10~ 


14 


ROSAT 

J. t KJ kJ /l J. 


600296N00 

KJKJKJ Z; VJL/ 


1992- 


-07-25 


PSPC 


2.5 


8.7 


x 


10~ 


14 


ROSAT 

1 1 KJ kJJi. J- 


600297N00 

KJKJKJ >J 1 ±y KJKJ 


1992- 


■07-25 


PSPC 


2 5 

Zi . W 


1.2 


x 


10~ 


13 


EOS AT 


600358N00 


1992- 


-07-26 


PSPC 


2.9 


4.1 


x 


10~ 


14 


ROSAT 

J. t yy kJ Ix. _L 


600303N00 


1992- 


■07-27 


PSPC 


2.4 


9.2 


x 


10~ 


14 


ROIAT 

1 1 KJ kJ J± J. 


600321 NOD 


1 992- 


07-29 

KJ t ZjC/ 


PSPC 


9 8 

Zi .O 


6 Q 


x 


1 0~ 


14 


ROIAT 

1 1 KJ kJ j1 ± 


600336N00 

KJKJ 


1 992- 


-07-29 

KJ t ZjC/ 


PSPC 


9 8 

Zi .O 


8 8 


x 


1 0~ 

-L KJ 


14 


ROSAT 

1 1 KJ kJ JrL -L 


600356N00 


1992- 


-08-03 


PSPC 


2.0 


2.5 


x 


IQ- 


14 


ROSAT 

J. t KJ kJ-ii. J- 


600338N00 

WVJVJOOOl'l KJKJ 


1992- 


■08-03 


PSPC 


2.2 


5.6 


x 


IQ- 


14 


EOS AT 


600361 N00 


1992- 


■08-05 


PSPC 


2.7 


4.8 


x 


IQ- 


14 


EOS AT 

l ( yy kJ-ii. -L 


600342N00 

WWWWTiiJ-L 1 WW 


1992- 


■08-05 


PSPC 


3.4 


6.8 


x 


IQ- 


14 


EOS AT 

i ( y y <j i i j 


600300N00 

UUU'J'JUl 1 WW 


1992- 


■08-05 


PSPC 


2.5 


8.6 


x 


IQ- 


14 


ROSAT 

1 1 KJ kJ jTL J_ 


600323N00 


1992- 


-08-05 


PSPC 


2.6 


1.5 


x 


IQ- 


13 


ROSAT 

J. t KJ J/1 J. 


600298N00 


1992- 


■08-05 


PSPC 


2.6 


2.2 


x 


IQ- 


13 


ROIAT 

1 1 KJ kJ J± J. 


600362N00 

KJKJKJ OKJZj±y KJKJ 


1 992- 


08-06 

UO KJKJ 


PSPC 


2.4 


5 7 


x 


1 0~ 


14 


ROIAT 

1 1 KJ kJ j1 ± 


60031 QN00 

KJKJKJ O .Lis ±y KJKJ 


1 992- 


-08-06 

KJO KJKJ 


PSPC 


1 8 


5 8 


x 


1 0~ 


14 


£?DQ A T 

liKJOJil 


DUUo4UlN uu 


1 QQ9 


Uo-UD 




z.D 


D.z, 


X 


1U 


14 


EOS AT 


600322N00 


1992- 


■08-06 


PSPC 


2.6 


1.0 


X 


10" 


13 


EOS AT 


600318N00 


1992- 


■08-06 


PSPC 


1.4 


1.2 


X 


10" 


13 


EOSAT 


600359N00 


1992- 


•08-06 


PSPC 


2.6 


1.7 


X 


io- 


13 


EOSAT 


600318a01 


1992- 


■12-21 


PSPC 


1.3 


1.2 


X 


10" 


13 


EOSAT 


600316N00 


1992- 


-12-21 


PSPC 


2.7 


1.4 


X 


10" 


13 


EOSAT 


600357N00 


1992- 


■12-30 


PSPC 


2.7 


4.1 


X 


10" 


14 


EOSAT 


600302N00 


1992- 


-12-30 


PSPC 


2.6 


1.3 


X 


10" 


13 


EOSAT 


600337N00 


1992- 


■12-30 


PSPC 


2.1 


1.5 


X 


10" 


13 


EOSAT 


600299a01 


1993- 


■01-01 


PSPC 


1.4 


6.2 


X 


10" 


14 


EOSAT 


600244N00 


1993- 


■01-02 


PSPC 


35.9 


1.8 


X 


10" 


14 


EOSAT 


600343N00 


1993- 


■01-04 


PSPC 


2.7 


7.5 


X 


io- 


14 


EOSAT 


600320a01 


1993- 


■07-01 


PSPC 


2.7 


6.3 


X 


10" 


14 


Chandra 


2046 


2000- 


■11-05 


ACIS-S 


14.8 


1.6 


X 


10" 


14 


Chandra 


2047 


2001- 


-03-06 


ACIS-S 


14.8 


2.1 


X 


10" 


14 


Chandra 


2048 


2001- 


■07-03 


ACIS-S 


14.0 


1.4 


X 


10" 


14 
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Table 3— 


-Continued 






Telescope 


ObsID 


Date 


Instrument 


rp a 
± exp 


Flux Upper Limit b 


XMM-Newton 


0402560101 


2006-06-28 


pn 


AAA 


4.2 x 10~ 15 


Swift 


00035337001 


2007-06-01 


XRT 


9.2 


9.0 x 10~ 15 


XMM-Newton 


0505760101 


2007-07-24 


pn 


57.0 


3.2 x 10~ 15 


XMM-Newton 


0511380101 


2008-01-02 


pn 


44.1 


1.5 x 10 12 


Swift 


00031919001 


2011-02-04 


XRT 


2.9 


3.3 x 10~ 14 


Swift 


00031919002 


2011-02-07 


XRT 


1.4 


3.3 x 10~ 14 


Swift 


00031919003 


2011-02-11 


XRT 


0.9 


3.3 x 10~ 14 


Swift 


00031919004 


2011-02-15 


XRT 


0.7 


3.3 x 10~ 14 



a Exposure time in kiloseconds. 

b Observed (i.e. absorbed) 0.2 - 10 keV 3a flux limits in units of erg cm~ 2 s _1 . For 
the 2008 XMM-Newton observation in which XMM J0038+40 was detected, we give the 
absorbed flux in bold derived from the double blackbody model. For the two remaining 
XMM-Newton observations, the upper limit was only calculated for the pn image. The 
images for the Swift data taken in 2011 were summed together to increase signal-to- noise, 
so the upper limits are for the combined data. While the upper limits are quoted for each 
individual ROSAT observation, the limit from the summed observations is 3.9 x 10~ 15 erg 
cm z s L . 
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